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Abstract 



The Hamiltonian limit of the corner transfer matrix (CTM) of a generalised free 
Fermion vertex system of finite size leads to a quantum spin Hamiltonian of the 
particular form: 

n N = - £ {* K< + i + + M< + <+i)) } 

n=l 

Diagonalisation may be achieved for all pairs of parameters (A, h) with the use of 
some new elliptic polynomials which extend the class of special polynomials known 
so far in the context of CTM. 



1 Introduction 



Recently we have studied the corner transfer matrix (CTM) of a free Fermion 8- vertex 
system at criticality This investigation has been part of a series of studies devoted 
to special cases in the parameter space of this model where an explicit solution can be 
found for the diagonalisation of such CTM's. The CTM is an interesting object in the 
statistical mechanics of two-dimensional lattice models which was introduced by Baxter 
in the seventies (see for a review and references to the original publications), and 
with which he was able to compute the mean magnetisation (or polarisation) in exactly 
solvable models, such as the eight- vertex model. Baxter has shown then that in the infinite 
limit of the two-dimensional lattice the spectrum of the "generator" of the CTM is an 
equidistant spectrum in some domain of the coupling constants. This result has been 
largely confirmed by numerous computations in vertex models of the RSOS type |3|. 

The pecularity of this spectrum has aroused interest in the last decade for another 
discovery in critical phenomena has a similar structure. The hypothesis of conformal 
invariance of critcal systems advocated by Belavin, Polyakov and Zamolodchikov [|J leads 
to the conclusion that the generators of row to row transfer matrices || as well as of CTM 
H of finite-size critical systems have also equidistant spectra. It is therefore natural 
to ask whether there exists a relation between these observations. In an attempt to 
clarify this question we have undertaken a systematic study of the CTM of the simplest 
soluble system, the free Fermion system (or its equivalent vertex model). We shall place 
ourselves usually at arbitrary temperature T and consider finite but large systems of size 
N. Therefore both limits T — > T c and N — > oo can be taken independently. 
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Let us summarise what has been obtained so far in the study of the CTM of the free 
fermion system. The most general free Fermion system depends on two parameters: the 
anisotropy parameter A, which measures temperature, and the reduced external magnetic 
field h 0. From the standard method of Lieb, Schultz and Mattis || it is known the CTM 
of a generalised free Fermion system may be diagonalised by diagonalising an associated 
matrix obtained from the eigenvalue equation of the generator of the CTM. This matix has 
the peculiar feature that its eigenf unctions are special polynomials defined by recursion 
relations. In reference || we have solved the simplest case: 

• A = 1 and h = 0, 

which is also a critial line. The polynomials obtained are particular cases of Meixner- 
Pollaczek polynomials. In reference |10| we have considered the case of 



• arbitrary A and h = 

corresponding to an Ising model and found two types of Carlitz polynomials of imaginary 
argument. Then several other cases are solved in f7|: 

• A = 1 and h arbitrary, 

where Pollaczek and Gottlieb polynomials are obtained, 

• A = h 2 , the disorder line UTTI, 



where Gottlieb polynomials are found. It is remarkable that the third class of Carlitz 



polynomials are also found in the CTM of the discrete Gaussian model [12]. Lately we 



have obtained also the generalised Pollaczek polynomials in a free Fermion vertex model 



with a line of defects Up to now the polynomials encountered are all orthogonal 

polynomials already known in the mathematical literature. Recently we have tackled 
other regimes in the parameter space (A, h) where the orthogonality property is not known. 
First we have studied in the polynomials associated with the critcal line 

A = 2h - 1 

and in this paper we shall deal with the general case 

\>h 2 

where both A and h are not restricted to satisfy any equation. The remaining regime of 
the parameter space, A < h 2 , where additional mathematical complexities arise, will be 
treated in a forthcoming publication. 

As the reader may have noticed the CTM of free Fermion systems is a simple device 
of statistical mechanics which introduces the special polynomials. The spectrum of a 
CTM of size N is "essentially" given by the zeros of a polynomial of size N, when N is 
sufficiently large. Hence, as one may guess, the distribution of zeros tends to be a uniform 
distribution as N —>■ oo. This fact confirms the findings of Baxter and verifies naturally 
the predictions of conformal invariance in critical systems. 

The new feature in this paper is the appearance of new "elliptic" polynomials which 
generalise the only known types of elliptic polynomials discovered by Carlitz fL4| . Some 
properties of these polynomials including the asymptotic distribution of their zeros, hence 
the eigenvalue spectrum of the system, will be given. 

Physically the "time" generator of such free Fermion 8-vertex problems is simply the 



anisotropic XY quantum spin chain in a magnetic field h. The counterpart of this using 
the CTM approach is the following generator for a finite chain of N sites: 

N-l 

Lo=J2 M«+i + A*Ku) + h(2n - 1)<} + h(N - (1) 

71=1 

where A is the anisotropy parameter describing essentially the temperature. Note the 
linear increase of the coupling strength along the chain. The problem of the simple chain, 
i. e. without that specific linear increase of the couplings, has been fully solved long ago 
| T5| . But the generator Lq has not yet, to our knowledge, been explicitly solved for general 
values of A and h. 

The paper is organised as follows. In section 2 we outline the method employed 
for calculating analytically the eigenvalues of the generator L which is based on the 
introduction of generating functions for the components of the eigenvectors of Lq. The 
recursion relations derived for the components of the eigenvectors are then equivalent to 
a set of coupled first-order differential equations which are solved formally in section 3. 
In the solution we encounter an elliptic integral. We need to appropriately parametrize 
and invert this elliptic integral which is done in section 4. This leads us to expressions for 
the generating functions in section 5 which are used in section 6 to derive explicitly the 
components of the eigenvectors as elliptic polynomials from Cauchy's theorem. In section 
7 finally we obtain an integral representation for the components of the eigenvectors 
that can be used to calculate asymptotically for large system size iV the eigenvalues of 
Lq. Section 8 summarizes our findings and gives an outlook on open problems. In the 
Appendix we give some details for certain limiting cases of the analysis of the main text. 
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2 Method for diagonalising Lq 

The standard method we adopt here is that of Lieb, Schultz and Mattis |§ which consists 
in rewriting Lq in terms of Fermion operators. Then one is left with the diagonalisation 
of two non-commuting matrices (A — B) and (A + B) in the language of Lieb, Schultz and 
Mattis. If we denote the components of the eigenvectors ip and by if) = (if)±, . . . , V'jv) 
and = (0i, . . . , 4>n) we have two recursion relations coupling the components if> n and n 

(n - l)V> n -i + n\ip n+1 - h(2n - l)ip n = e<j) n (2) 
A(n - l)0 n _i + n<p ri+l - h(2n - 1)0„ = eip n (3) 

For the end components n = N of the finite chain we have only 

(N - l)^_i - h(N - l)ifj N = e<p N (4) 
X(N - 1)0^! - h(N - l)<f) N = eifj N (5) 

Using the recursion relations @) and (f|) the equations for the end components (£|) and 
(|) are equivalent to 

0v+i = h(f> N (6) 
\if) N+ i = hip N (7) 

which resemble periodic boundary conditions. Note that the components (pN+i and ipN+i 
are only defined through these equations. 

The method employed is simple in principle: it consists to finding an expression for 
ip n and n , then on substituting into eqs. ([|) and (^) one obtains the eigenvalues e. 
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However, as we are only interested in large values of N an asymptotic expression will only 
be necessary. For the so far known cases pj the ip n and <p n are expressible in terms of 
known special polynomials. We expect thus that for general values of A and h they are also 
polynomials of the elliptic type which are seen in the case of the doubled Ising model |10| . 



In special cases where A and h were chosen to satisfy particular conditions it was previously 
possible to identify the recursion relations (§) and after some trivial transformations, 
to the recursion relations for some known polynomials. Here as in another case discussed 
recently |13] we introduce the generating functions as formal power series in a parameter 



t with the components of the eigenvectors as coefficients 

oo 

m = E *""Vn (8) 
n=l 

oo 

= 53 r-Vn (9) 

n=l 

The recursion relations (|2|) and (|3|) are then equivalent to a set of coupled first-order 
differential equations for the generating functions 

(t 2 + A - 2ht)ip' + (t- h)ijj = E<p (10) 
(At 2 + 1 - 2ht)4>' + (At - h)(p = e^j (11) 

From the explicit solutions ip(t) and 0(t) of these differential equations with the given 
initial conditions one may then extract <p n and ip n by using the Cauchy theorem. Then for 
large N the boundary conditions (^) and (|7D determine the spectrum of L as in previous 
works of this series Jl|, 0]. Since the presence of two parameters A and h complicates 
the mathematical working considerably we shall proceed step by step in presenting the 



solution. The central problem at hand is simply the parametrization by elliptic functions 
of the solution of the coupled set of differential equations (pTOj ) and fllTD - 



3 Formal solution of the differential equations — gen- 
erating functions 

Before writing down the formal solution of the differential equations (|1(]) and ([11]), we 
note that these equations as well as the recursion relations for ip n and </>„, eqs. (^) and 
remain globally invariant under the combined transformations 

A — > A -1 , h — > \^ 1 h and e — > A _1 e 

followed by 

4> -> and ^ ^ (4>n -> 0n and n -> Vn) 
This property allows us to restrict our study to the domain S defined by 

< A < 1 and < h < 1. 



Using this symmetry relation integrating factors can be found |16| and as in [|IJ the general 



solutions of the differential equations take the form of Meixner's generating functions [pjj 
(up to constants) 

if>(t) oc fit) exp(£w(t; A, h)) (12) 
4>{t) oc g(t) exp(ew(t; A, h)) (13) 



where 



f(t) = [t 2 + \-2ht)—2, g(t) = (At 2 + 1 - 2ht)—2 (14) 



and w = w(t; A, h) is the two-parameter elliptic integral 

w(t;X,h)= / dxf(x)g(x) (15) 
The lower bound of the integral will be chosen as to agree with known results. Note that 



for h = 0, one recovers the case of the doubled Ising model studied in [10|. Generating 



functions of the type flT2]) and (|13|) are called generating functions of the Meixner type 



by Chihara ||18|| . Thus it appears that (p^) and (p~3| ) represent perhaps the most general 
case known so far. The main problem here is the inversion of the two-parameter elliptic 
integral w(t), eq. QToD, and then the representation of the ip n and <fi n as contour integrals. 

4 Parametrization and inversion of the elliptic inte- 
gral w(t) 

An essential step in the explicit solution consists in an appropriate parametrization of the 



elliptic integral w(t). We shall follow here the procedure given by Greenhill |T9|. Both 
polynomials N(t) = t 2 + A — 2ht and D(t) = At 2 + 1 — 2ht appearing in w(t) have the 
same discriminant A = h 2 — A. The zeros of N(t) in the parameter range < A < 1 and 
< h < 1, the region S, are 

y = h + VK, 6 = h-VK (16) 

whereas D{t) has zeros at 

<*=l> (17) 

The square 5* is divided into 2 regions by the disorder line J?| A = h 2 (see Fig. 1), 
separating oscillating from monotonous behaviour of the correlation functions JTl . 



In the region Si, A > h 2 , the zeros are pairwise complex conjugate, whereas in the 
region S 2 , A < h 2 , the zeros are all real and ordered according to— oo<<5<7</3< 
a < oo. 

The inversion of the elliptic integral in the region S2 is more involved and we shall 
defer its study to a sequel publication. In the remainder of this paper we shall be only 
concerned with the region Si. 

The region Si is limited by three boundaries which contain known results @: 



• h = The doubled Ising model, viewed as free Fermion eight- vertex model |L0[ : ip n 
and (f) n are Carlitz elliptic polynomials of imaginary arguments. 

• A = 1 Isotropic case in the presence of a magnetic field h. The solutions ip n are 
given in terms of Meixner Pollaczek polynomials. 

• A = h 2 The disorder line where the solutions ip n are expressed in terms of Gottlieb 
(Meixner polynomials of the first kind) polynomials. 

Note that on these three boundaries the three types of polynomials are all orthogonal 
polynomials, their recursion relations are always reducible to tridiagonal form. The ip n 
and 4> n studied here provide an interpolation between the three classes of polynomials and 
are presumably also orthogonal polynomials. 

The elliptic integral w(t) may be inverted in a standard way: one may give t as a 



function of w following the example given in the book of Greenhill [19|, §70. In the 
following we outline Greenhill's procedure. 
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Through the change of variables 

N(t) t 2 + X — 2ht 
V{t) ~W)~ W + l- 2ht (18) 

the elliptic integral is transformed from its Jacobian form w(t) into the WeierstraB form 

. . 1 [V dx , , 

w{y) = ^==^ m (19) 

VA - h z J yo ^Ax(y 1 -x)(x- y 2 ) 

where yo = y(to) and y± and y 2 are the maximum and minimum values of y(t), respectively, 
at the points t\ and t 2 related by 

^ 2 (l-yi, 2 ) 2 = (A-y 1)2 )(l-Ay 1)2 ), y, = y, 1 (20) 

^2-^^^,2 + 1 = 0, h = t 2 l (21) 

and 

V2 = 2{x l _ h2) {(1 + A 2 - 2h 2 ) - (1 - A)^/(l + A)2-4^} < yi (22) 

Then from the standard WeierstraB form one may solve the inversion problem in three 
different ways 



w(y) = J^sn- 



lyl- 


-y 


yi - 




Iv- 


1/2 



! ' 2 rn : I . - —.i< I (23) 



A-/i 2 \V yi — S/2 



^ :dn^( 



V A-/i 2 

in terms of the Jacobian elliptic functions sn(v,K), cu(v,k) and dn(i>, k) of argument v 
and modulus k. In particular the modulus k is given by k 2 = 1 — y%, thus 1/2 is the 
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complementary modulus of k. Conversely for later use one can extract from (^3|) the 
quantities 

2/1- V = (Vi - V2)sn 2 (wq, «) 

y -yz = (yi - y 2 )cn 2 (wq,K) (24) 
y = — dn 2 (wq, k) 

2/2 



with q = \J—^-- To obtain i as a function of w, i. e. to invert the elliptic integral fll5|) , 
we use equation (102) of Greenhill |19| 



Vl ~ y = At 2 + 1 - 2ht (25) 

y-m - (1 - A ^ 2)(t - t2)2 (26) 

y y A/ 2 + 1 - 2/it v ; 



By dividing out these equations and using the first two of eqs. ( P^ ) as well as y\ ■ y 2 = 1 
and t\ ■ t 2 — 1, we obtain the following form for t = t(w) 



, , *i - y^cn(qw, k) + t 2 Jy 2 {l - Xy 2 )sn(qw, k) 

t(w) = — = / (27) 

VA - y 2 cn(qw, k) + Jy 2 (l - Xy 2 )sn(qw, k) 



Thus in the case of two parameters A and h the inversion of eq. ([15]) yields a rational 
function of the Jacobi elliptic functions sn and cn. However in order to agree with the 
known results on the boundary line h = 0, where we have a canonical elliptic integral of 
Legendre form and hence t is simply proportional to a Jacobi elliptic sn(v , k) function, 
as was always the case in previous studies J/J on special lines where one had therefore 
only one parameter, we transform to imaginary argument and complementary modulus 

2/2 (ft 2 + 2/1 = 1) with the transformation formulae (Jacobi's imaginary transformation) 

sn(iqw, k') 1 

sn(qw, k) = -i — -, cn(qw, K) = — (28) 

cn{iqw, k') cn{tqw, k') 
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Thereby we obtain 



1 - it ly /y^n.{iqw,y 2 ) 



At h = one may directly invert the elliptic integral of Jacobian form 

i doc 

w(t) 

to obtain 



o ^(Aa^ + lX^ + A) 



. . /— sn(iy.A') /- - 

t( w ) = VA • A' = Vl-A 

cn(ty, A') 

By the transformation to imaginary argument and complementary modulus this yields 

t = —ivXsn(iw, A) (30) 



Of course we want that (p9|) reduces to ( p0|) in the limit fa — > 0. To achieve this we shift 
the argument in (^) by iT', the complete elliptic integral of the complimentary modulus, 
i. e. the quarter period in the imaginary direction of the Jacobian elliptic functions: 
w = —u + Wo with qwo = K'tyz), which transforms the function sn according to 

Through this last transformation we arrive at 

h - iJy~2~sn(iqu, y 2 ) 



, , _ . , (31) 

Basically this result means that we have to choose the lower integration bound in ([15]) in 
such a way that ( |3~TD holds. For h — > (i. e. y 2 — > A and t 2 — > 0) now ( pi]) agrees with 
(BO)- From now on we shall use this condition at h —>■ as reference which will also fix 



the constants in the generating functions (|T^D and fllci) 
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5 Expressions for the generating functions 



From (|26|) we have for the prefactor of the generating function (j){t) to be 



VXt 2 + l-2ht 



Xy 2 



{t - h) 



V V-V2 

Using the second of eqs. fl24f) and performing both transformations described above con- 
secutively, i. e. Jacobi's imaginary transformation and the transformation according to 
the shifted argument, we obtain together with fl3~i~D to evaluate (i — t 2 ) 

1 



Vl t-ihv / y2sn(iqu,y 2 ) 



VXt 2 + l-2ht \ 1 - Xy 2 (1 - tl)dn(iqu, y 2 ) 



For h — > this expression has the correct limit [ItJ 

lim . = — — — 

h->oo VAt 2 + l- 2ht dn(iqu, X) 

Hence we ontain the generating function of the n as 



VXt 2 + l-2ht 



1 — 2/1 exp(g 1 eK'(y 2 )) 1 - it 2y /y2sn(iqu, y 2 ) 



\\l-Xy 2 

Similarily using the first of eqs. 



'1-t 2 ) 



dn(iqu, y 2 ) 



(32) 



and eq. (j25|) we obtain after the necessary trans- 
formations the generating function for the other set of components of the eigenvectors 



"0n 



Vt 2 + A - 2ht 



y\ exp(g 1 eK'(y 2 )) 1 - it 2 Vy^sn(iqu, y 2 ) 



\(l-Xy 2 ) (1 - 1|) 



y 2 cn(iqu, y 2 ) 



(33) 
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To simplify the notation we shall set 

AT 

and 



exp(g eK (y 2 )) 



\0- Ay 2 ) 
f = iy / y 2 'sn(iqu,y 2 



so that 
<P(t) = 



M l ~ U -e^ 1 ^, and m 



l-t|) dn(iqu : y 2 ) 



l-tV, 



y 2 cn(iqu, y 2 ) 



(34) 



(35) 



iq 1 e(iqu) 



(36) 



6 The elliptic polynomials ip n and 



n 



From the generating functions (|T2|) and (12) we may compute the coefficients ip n and 
which are of course functions of the eigenvalue e with the help of Cauchy's formula 

= —flf>(t)t- n - 1 dt 



Substituting vp(t) by its expression ( p6|) and using the variable £ with 



we obtain 



t 



J\f 



dt 



:i-^) 2 

(i-t 2 £) n 



2m 7 y/y^ca(z,y 2 ) (h - n+1 
where z = iqu and x = q~ l e, the scaled eigenvalue. 

We can express tf) n in terms of the usual Carlitz polynomials (cf. |T 

e ixz _ ~ f e y p;{x) 



(37) 



(3? 
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where 

P* (t \ = j C p for P even 
p { ' \ D; for p odd 

are the Carlitz polynomials of even and odd order, respectively, associated with the Ja- 

cobian elliptic functions cn and dn. Then ip n takes the form of a formally infinite series 

in these polynomials 

Mx) = _^ I . W . (39) 



The condition i = is equivalent to ( = t 2 — £ = 0. Thus the contour integral in 
may be directly evaluated in the complex ( plane 



(*a " C) P - ^ = 2ma n {p) 



where 



We arrive at 



*n(p) = E (-ir* fa)fm) (1 - t 2) m fe)- m+( ' (40) 

q+m=p \ " / \ / 



* w ^£- (p) (Jfv (41) 

However the ip n {x) are nevertheless polynomials of finite order n in the variable x. To see 
this we first evaluate formally for n = 0. We have 

^o(p) = 



and 

p 
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Defining z 2 by 
we can resum the series and find 



Kf e ixz 2 



which is essentially a constant. From ipo taken formally at arbitrary values of z we can 
compute all ip n for n > and thereby indeed show that ip n are polynomials of finite order 
n. For ipi we obtain e. g. the following formula 

Mx) = {t2-{i-4)-^Mx) 

which is a polynomial in x of first order after dividing out the exponential. 

More generally we can show that ip n is a polynomial of order n in x explicitly given 

by 

"' + ^(n-M^'W- <43) 
With the same procedure we show that <p n are polynomials in x of finite order n 

^-"S^^r (44) 



Here 



_ / d p for p even 

V P l x) I c * for p odd 
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and we have for n = 

J.XZ2 

0o(x)=A/ — r (45) 

dn(z 2 ,y 2 ) 

from which we again compute all polynomials (j) n {x) of order n > explicitly 



7 Real integral representation and asymptotic be- 
haviour 



In the same spirit as in reference |T0|| we derive now the real integral representation for 



^pn(x), using eq. (p7\) in the complex z-plane. As before in |T(| we choose the rectangular 
contour lm(z) = ±K'(y 2 ) and He(z) = ±K(y 2 ) surrounding the point z 2 which now is not 
necessarily at as it was the case in [jID| |. We obtain ip n {x) as a sum of two contributions 



of the form 



y/V^N [ K{V2) \ aX K'(y 2 )+i X v (h + i^msn(v,y 2 )) n 



i'n(x) = / < e 



2tt J-K(y 2 ){ (l+it 2y /y 2 ~sn(v,y 2 )) n+1 
(I - it 2 ^/y 2 ^sn[v,y 2 )) n + 1 ) 



KJ 



' N f K(K) SjxK'M-xv (t2y/m + idn(v, K )) 



2tt J-k(k) { ( v /^ + zt 2 dn(u,/«)) n + 1 

(V^-«t 2 dn(v, «))"+!/ V ; V ; 
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Note that in the second part of eq. ([47]) the elliptic modulus k — y 1 — y 2 appears, the 
complementary modulus of y 2 . As in [|H]] the second contribution in ( f4"T| ) dominates when 
n — > oo because the complex number 

t 2 ±i y /y 2 ~sn(v,y 2 ) 



1 ± it 2 ^sn(v,y 2 ) 



has a modulus smaller than 1 



t\ + y 2 sn 2 (u,7/ 2 ) 



< 1, 



\| 1 + ti?/2sn 2 (w,7/2) 

since ?/2sn 2 (t> , y 2 )) < 1 is always fulfilled for general values of v. On the contrary the 
modulus of the complex number 

t 2y /y 2 '±idn(v, k) 



y/y^ ± it 2 dn(v, k) 
may be larger than 1, since y 2 < dn(t>, k) < 1, i. e. 



t\y 2 + dn 2 (v,K) ( > 1 for < v < K/2 
\y 2 + t 2 2 dn 2 (v, k) 1 < 1 for AT/2 < v < K 

Hence for n->oowe only have to consider 



ixK' (k)—xv 



{t 2 ^ + idn(v, 



+e 



2tt J-k(k) { 

-ixK'( K )+xv (t 2 y/yi-idn(v,K)) 



y 2 + it 2 dn(v, n)) n+1 
cn(v, n)dv 



(4£ 



(^/m-ihdniv,^)^ 1 

Using that dn(v, k) and cn(v , k) are both even functions of the variable v, we derive an 
alternative form for the asymptotic of ip n (x) 

2kM 
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{cos(xK' (K))I n (x) — sm(xK' \n))Z' n (x)} 



(49) 
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where T n {x) and 1' n {x) are the following integrals 

i-k(k) ^Jy^cosmp + t2dn(t> , k) sinmp 

y 2 + tldn 2 (v, k) 
k (k) ^/y^sinrup — t 2 dn(t>, k) cos nip 



1 n \x 



- L 

JO 



1/2 + £ 2 dn (v, k) 



cosh(a;v)cn(v, n)p n {v )dv (50) 
cosh (rev )cn(v, n)p n (v)dv (51) 



with p{v ) and <p(v) defined by 



p(v)e 



iip(v) 



hy/yl + idn(v, k) 
v /y 2 " + it 2 dn(v, k) 



(52) 



We observe that p(v ) is a decreasing function of v, aquiring values from \j "^Js" at v = 



to 



V2+4 
t%V2 + l 



at v = K(k). Since p(\K) = 1 the integration intervall for large n — > oo 



is practically limited to [0,1^/2] only and we can approximate p by its maximum value 



P 



U2+i. 



t. Moreover in this interval, cosh(a;v)cn(t;, k) = 1 and peaks at a point near 



v = K. This means that T n (x) = X n and T' n (x) = T' n are constants independent of x 



%\)2 + 1 i f\ K ( K ) ^fy~2 cos rup + t 2 dn(v, k) sinn<p 



y 2 + tldn 2 (v, k 



XL 



t\y 2 + 1 \ t\ K ( K ) ^Jy 2 s'm.wp) — t 2 dn(t>, k) cosnip 



dv 



dv 



(53) 
(54) 



J 2/2 + ^1 I Jo y 2 + 4dn 2 (v,n) 

Within these approximations one may estimate the zeros of i[) n (x) m the asymptotic limit 
n — > oo. They are given by the equation 



tfm(xK'{n)) 



or, for each n, the zeros x np are labeled by p 



T 

J -n 

u 



tan 9 % 



(55) 



x np = 9 n +pir 



(56) 
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The eigenvalue problem set in eqs. (|2] - [7|) is equivalent to 



tan(x Np K'(K)) = — = tan 6^ (57) 

a± n+1 n.L N 



leading to 



V V2 K'(k) 



The spacing between the ejy p is thus constant and there is a shift for each N given by 6^ 
defined by eq. (pwj). 

A similar analysis for <p n (x) may be given. We only quote the results 



+e 



2tt J-K( y2 ) { (1 + it 2y /y^sn(v, y 2 )) n+l 

+ mAf /•*(«) [ cixK 'M- X y (t2^ + idn(v,K)) n 
2n J-K( K )\ ( v ^ + it 2 dn(t»,/t))«+ 1 



-ixK'{K)+xv {t 2y /yj-ifoi{v,K)y 
fyi - it 2 dn(v, 



lxK w+xv v^v^ v ; a/ , ^ / 5g x 



Again only the second integral dominates as n —>■ oo for the same reason as before in the 
expression for ip n . This last part may be recast into the following form which incidently 
is proportional to i 



71 



with 



-{cos{xK'(K))I n (x) - sin(xK'(K))l' n (x)} (60) 



Z n {x) = [ y^cosmp + t 2 dn{v, k) smnip ^^j^^ K)p n {v)dv (61) 

Jo y 2 + t 2 dn [v, k) 

T n {x) = f Vy^ smni P J^ V,K °° Snip smh(xv)sn(v,K) p n (v)dv (62) 

Jo y 2 + t 2 dn (v, k) 



21 



where p and ip are again given by ([52]). For n — > oo the behaviour of p limits again 
the integration to the interval [0, K/2] and p can be approximated by its maximum value 
p = J ~^^r - However here the product sinh(xf )sn(w, k) behaves as xv 2 in [0, K/2] instead 
of beeing nearly a constant of order unity. Thus we obtain the asymptotic behaviour for 

(f>n(x) 



4>n(x) 



2%kN 



71 



{x cos(xK \n))X n — x sm(xK' (K,))X' n } 



(63) 



where now 



T 



tlV2 



\ V2 + t 



tlV2 



2 \[yi cos n( P + t 2 da(v, k) sin ncp 

V — „ . o , ; aV 



y 2 + *§dn (v, k) 



2 \/lE sin nip — t 2 dn(v, n) cos nip 
v — = , o . : dv 



J V2 + % ) Jo y 2 + tldn 2 (v,K) 

are constants independent of x. The of 4> n (x) are asymptotically given by 

ta.n(x np K' (k)) = -=^- = tan# n 

The boundary conditions @ and (^) of the CTM problem yield the eigenvalues tN P 

Xtv+i — hX N 



(64) 



tan(g l e^ p K' {k)) = tajo.0 



N 



V — hV 
-L N+ i n-L N 



(65) 



or 



K'(k 



■{9 N +prr} 



(66) 



We observe again that the levels e^ P are equidistant in this limit of large N, but there is 
a translation by an amount k ?i k \ Qn for each polynomial of order N. 

Eqs. fl58|) and ( |66"D are the main results of this paper, again confirming, for the simplest 
vertex model, the results of Baxter B in an explicit calculation for a finite system. 
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A last remark concerns the level spacing. From (|58| ) and (^) we have 

Ae = e N , p+1 - e N , p = = ^ (67) 

The modulus used here is 2/2 of eq. (^). For /j = 0we have 2/2 = A which is thus different 
from the modulus parametrization used in which was A -1 . Moreover the normalization 
of the energy levels is also different due to the fact that the Carlitz polynomials are directly 
used in ||10|| . The Carlitz polynomials are normalized so that the first one is always 1 or 



x depending upon the parity. Here our ipo{ x ) is n °t 1 but contains an x dependence 
according to eq. (f4~2|) which may be divided out later. 

At h = where 2/2 = A and t 2 = there is a decoupling in eqs. ([47]) and (|59|), 



respectively, and we recover the results of |L0 



As in [|IJ we have checked the level spacing ( |6^j ) numerically with standard methods 
(cf. |2l|) of diagonalising the pentadiagonal matrix which is equivalent to the recursion 
relations (0) and (Q) if either ip n or <p n is eliminated. We observed equidistant level spacing 
to rather high accuracy already for very moderate system sizes of the order of iV = 20. 

8 Summary and conclusion 

The generator of the CTM of a generalized free Fermion vertex system of finite size is a 
quantum spin chain Hamiltonian with particular interactions which increase linearly along 
the chain. We have presented the analytical diagonalisation of this particular quantum 
spin chain in the asymptotic regime of large system size N for arbitrary values of the 
parameters, the anisotropy A and the magnetic field h, in the region where A > h 2 . 
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Let us briefly summarise the methods applied to accomplish our goal and restate our 
main result for easy reference. The asymptotic diagonalisation has been achieved through 
the explicit construction of a new class of elliptic polynomials which are the components 
of the eigenvectors of the problem. In this construction an elliptic parametrization of the 
generating functions of the polynomials has been used which is based on the treatment of 
a two-parameter elliptic integral. The asymptotic evaluation of an integral representation 
of these polynomials yields the eigenvalues, given by eqs. ( |58|) and (|66D , respectively, which 
are equidistantly spaced with spacing 

A ^ 

Ae = ^ 

the modulus yi of the complete elliptic integral K(y 2 ) being related to the generating 
functions of the eigenvectors and given explicitly in (|22|). This equidistant spacing is 
the main result of the present work extending the findings of previous studies [Q, [7|, |], 
|7| to general values of the parameters and thereby confirming once again the general 
expectation 0. 

We have not touched the issue of the orthogonality of the polynomials ip n (x) and 4> n (x) 
which may be called associate Carlitz polynomials. Since the three limiting cases are 
orthogonal polynomials, it is natural to expect that the "associated Carlitz polynomials" 
remain orthogonal. Presumably the proof is based on the continuous fraction expansion 
of some elliptic functions interpolating between the Jacobian elliptic functions cn(x, k) 
and dn(x, k) [TJJ]. We have not succeded in proving this yet. 
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Appendix 



As we have seen the main problem has been the inversion of the elliptic integral of [Lq . 
The solution obtained was given in eq. (|3ll) which upon replacing t 2 by its expression 
reads 

f _ h(l - 2/2) - iQ- ~ Ay 2 )y/y2sn(igM,|/2) 
(1 - Xy 2 ) - ih(l - y 2 )^sn(iqu,y 2 ) 

In this appendix we check the limit A — > 1 against a direct calculation. Let us assume 
A = 1 — e, then y 2 = 1 — yj^f or using h = cos as in ref. y 2 = 1 — ej sin^. Then to 
first order in e (1 — Xy 2 ) = e(l + 1/ sin6>), and we have 

lim sn(iqu, y 2 ) = tanh(ra sin 6) 

A — > 1 

and 

cos 9 + (1 + sin 9) tan(u sin 9) 



lim t 



A-»l (1 + sin^) + cos^tan(wsin^) 
which upon inversion gives u as a function of t 

1 /(l + sin0-icos0) +i(cos0- (1 + sin0)t) > 
U ^ ' ~ 2i sin 9 U \ (1 + sin 9 - t cos 0) - i(cos - (1 + sin 9)t) 
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or 

1 ( t-e~ ie \ + 7T/2 
M( } ~ 2ism9 \t-e i6 ) 2sinfl 

which agrees with the u(t) directly computed from the integral ([15]) with A = 1. 
Next we check the limit A — > h 2 . There we have 

X-h 2 

and 



lim i^/y^sn(iqu, y 2 ) = - ^ sinh((l - h 2 )u) 
2/2—0 — h A ) 

Now to obtain the correct limit we must impose a shift 

u(l - h 2 ) = x(l - h 2 ) - l -K\y 2 ) 

Then as y 2 — > 



lim = x(l — h ) InU/yo) = x(l - h ) - In . , 

?/2-o v ; 2 v/ y ; v ; \ V X-h 2 

leads to 



Finally for A — > /i 2 



sinh((l-/,»^-i^=§e-^ 2 ) 



lim iJy2~sn(iqu,y 2 ) = e x{1 h } 

2/2—0 



which is want one obtains by direct integration. 

The correctness of the two limits implies that the limiting generating functions are 
generating functions for Meixner polynomials of first and second kind according to the 
classification of Chihara |TB[ or the Gottlieb and Meixner Pollaczek polynomials according 
to an independent classification. 
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